Fourth-order Runge-Kutta Method Research Articles

Effect of Interstage Torsional Damping on Nonlinear Dynamic Characteristics of Two-Stage Planetary Gear System

Abstract This paper explores the influence of interstage torsional damping on the nonlinear dynamics of a two-stage planetary gear transmission system. A comprehensive nonlinear dynamic model is developed, incorporating interstage torsional damping, stiffness, time-varying meshing parameters, damping, and tooth side clearance. The nonlinear equations are derived and solved using the fourth-order Runge–Kutta method. The investigation reveals that the changes in torsional damping between stages significantly affect the nonlinear dynamic characteristics of the system, especially under load fluctuations. Bifurcation characteristics are analyzed using phase diagrams, Poincaré diagrams, time history diagrams, and frequency–domain spectrograms. The simulation results demonstrate that heightened torsional damping between stages distinctly affects the motion states across excitation frequencies. Notably, it mitigates the unstable motion state caused by high-frequency excitation when tooth side clearance is minimal.

Study of the Six-Compartment Nonlinear COVID-19 Model with the Homotopy Perturbation Method

The current study aims to utilize the homotopy perturbation method (HPM) to solve nonlinear dynamical models, with a particular focus on models related to predicting and controlling pandemics, such as the SIR model. Specifically, we apply this method to solve a six-compartment model for the novel coronavirus (COVID-19), which includes susceptible, exposed, asymptomatic infected, symptomatic infected, and recovered individuals, and the concentration of COVID-19 in the environment is indicated by S(t), E(t), A(t), I(t), R(t), and B(t), respectively. We present the series solution of this model by varying the controlling parameters and representing them graphically. Additionally, we verify the accuracy of the series solution (up to the (n−1)th-degree polynomial) that satisfies both the initial conditions and the model, with all coefficients correct at 18 decimal places. Furthermore, we have compared our results with the Runge–Kutta fourth-order method. Based on our findings, we conclude that the homotopy perturbation method is a promising approach to solve nonlinear dynamical models, particularly those associated with pandemics. This method provides valuable insight into how the control of various parameters can affect the model. We suggest that future studies can expand on our work by exploring additional models and assessing the applicability of other analytical methods.

Generalized Shannon entropy sparse wavelet packet transform for fault detection of traction motor bearings in high-speed trains

An effective structural health monitoring method of traction motor bearings is a powerful guarantee for the safety operation of high-speed trains. However, it is exceptionally difficult to detect bearing fault characteristics from the vibration signals of traction motor bearings operating at high rotational speeds. In this scenario, a generalized Shannon entropy sparse wavelet packet transform (GSWPT) for fault detection of motor bearings is proposed in this paper. Firstly, a generalized Shannon entropy sparse regularization method is proposed to obtain sparse wavelet reconstruction coefficients by extending the definition of the Shannon information entropy, and the non-convex sparse regularization function is minimized by synergistic swarm optimization algorithm. Then, the wavelet node coefficients are weighted according to the second-order cyclostationarity index of the wavelet packet node to further enhance the sparsity of the reconstructed signal. Moreover, the optimal decomposition level of GSWPT is adaptively selected by the maximum sparsity and cyclostationarity criterion. Particularly, in order to verify the bearing fault detection performance of GSWPT in practical engineering, a bearing fault dynamic model of traction motor in high-speed train was established based on Hertz contact theory and the fourth-order Runge-Kutta method to obtain simulated data under strong Gaussian white noise, and a corresponding test platform was constructed to collect experimental data under different operating conditions. Finally, the applications on the simulated and experimental signals of traction motor bearings in high-speed trains demonstrate that GSWPT significantly outperforms the conventional wavelet packet transform, dual-tree complex wavelet packet transform, blind deconvolution, modal decomposition, and Infogram methods to some extent for fault detection.

Nonlinear vibrations, bifurcations and chaos of piezoelectric composite lattice sandwich plate with four simply supported edges

This paper is devoted to investigate the nonlinear dynamics, bifurcation and chaos of the piezoelectric composite lattice sandwich plate with 1:3 internal resonance. Galerkin method is used to discretize the nonlinear partial differential governing equations of motion into the ordinary differential equations. The multi-scale method is used to obtain the modulation equations of the piezoelectric composite lattice sandwich plate in the polar coordinates and Cartesian coordinates in the primary and parametric resonances, respectively. Using Newton-Raphson method, the bifurcation diagrams of the steady-state equilibrium solution are obtained for the modulation equation with the change of the system parameters. The stability of the steady-state equilibrium solution is analyzed for the modulation equation. The existence of the static bifurcations, such as the saddle-node bifurcation, pitchfork bifurcation and Hopf dynamic bifurcations are found. The complex nonlinear jump phenomena caused by the existence of the multiple nonlinear steady-state solution regions are analyzed in detail. Fourth-order Runge-Kutta method is used to continue tracing the nonlinear periodic solution of the modulation equation in Cartesian coordinates. The time-history diagram, phase portrait and Poincaré map are obtained for the piezoelectric composite lattice sandwich plate in the nonlinear resonance. The dynamic process of the nonlinear periodic solution is described for the modulation equation from the limit cycle oscillation to the chaotic attractor in detail.

Design and optimization of a two-degrees-of-freedom single-sided vibro-impact nonlinear energy sink for transient vibration suppression of a thin plate

To rapidly suppress the transient vibration response of a thin plate, a two-degrees-of-freedom single-sided vibro-impact nonlinear energy sink (TSSVI NES) is proposed in this paper. Based on the Kirchhoff hypothesis, the coupled equations of motion between the plate and the TSSVI NES are established. The system is modeled by the Galerkin technique and numerically solved using the fourth order Runge-Kutta method. To obtain the optimal parameters of the TSSVI NES, a particle swarm optimization method is adopted. The system response with the TSSVI NES is compared with that of the cubic nonlinear energy sink (CNES) and the single-sided vibro-impact nonlinear energy sink (SSVI NES), showing that the vibration suppression ability of the TSSVI NES is better. The mechanism of rapid response suppression of the TSSVI NES is demonstrated by way of simulations. A parametric study is also conducted to determine the effects of damping, the coefficient of restitution, mass, stiffness, and installation position of TSSVI NES on the efficacy of the vibration control device. The numerical results demonstrate that the performance of the TSSVI NES is robust to the damping, coefficient of restitution, and installation position, but sensitive to mass difference and stiffness difference. Finally, the effect of the TSSVI NES is examined when the plate is excited by seismic loads, demonstrating a good response suppression performance.

Study on the steady state position of ocean towing cables

A study is conducted on the static equilibrium posture of the towing cable in water under a constant speed state in the towing system. Firstly, the underwater force state of the towing body is calculated and then the posture and tension of the towing cable are analyzed and calculated. As the cable is divided into two states which are above water and below water, a differential equation system is established as the boundary. The initial value of the equation can be transformed by the drag and gravity values of the towing body. The program solves the differential equations using the fourth-order Runge Kutta method and draws the attitude curve. Finally, it verifies the theoretical calculation results through experiments.

Numerical Simulation on the Two-Degree-of-Freedom Flow-Induced Vibration of a Submerged Floating Tunnel under Current

The submerged floating tunnel (SFT) is a novel form of transportation infrastructure for crossing deeper and wider seas. One of the primary challenges in designing SFTs is understanding their hydrodynamic response to complex environmental loads. In order to investigate the two-degree-of-freedom (2-DOF) flow-induced vibration (FIV) response of SFTs under current, a two-dimensional (2D) numerical model was developed using the Reynolds-averaged Navier–Stokes (RANS) method combined with the fourth-order Runge–Kutta method. The numerical results were validated by comparing them with the existing literature. The study then addressed the effects of coupled vibration and structural parameters, i.e., the mass ratio and natural frequency ratio, on the response and wake pattern of SFTs, numerically. The results indicated that coupled vibration had a significant impact on the SFT response at reduced velocities of Urwx ≥ 4.4. A decrease in mass ratio (m* < 1) notably amplified the 2-DOF vibration amplitudes of SFTs at Urwx ≥ 4.4, particularly for in-line vibration. Similarly, a decrease in natural frequency ratio (Rf < 1) significantly suppressed the in-line vibration of SFTs at Urwx ≥ 2.5. Therefore, for the design of SFTs, careful consideration should be given to the effect of mass ratio and natural frequency ratio on in-line vibration.

Two-dimensional Coupled Asymmetric van der Pol Oscillator

A reliable algorithm based on an adaptation of the standard differential transform method (DTM) is presented, which is the multi-step differential transform method (Ms-DTM) since it may be difficult to directly apply differential transform method (DTM) to obtain the series solutions for the present two- dimensional coupled asymmetric van der Pol oscillator. The solutions of a two dimensional coupled asymmetric van der Pol oscillator were obtained by Ms-DTM. Figurative comparisons between the Ms-DTM and the classical fourth order Runge-Kutta method (RK4) are given. The obtained results reveal that the proposed technique is a promising tool to solve the considered van der Pol oscillator and yield same information on the phase portrait confirming the stability of the system, effectively. It can be said that the considered approach can be easily extended to other nonlinear van der Pol oscillator systems and therefore is widely applicable in engineering and other sciences.

Impact of Silver Nanoparticles Infused in Blood in a Stenosed Artery under the Effect of Magnetic Field Imp. of Silver Nano. Inf. in Blood in a Sten. Art. Under the Eff. of Mag. Field

Nanotechnology and nanofluids offer promising avenues for enhancing our understanding and control of blood flow dynamics, providing novel solutions for addressing challenges in cardiovascular health and medical diagnostics. This paper investigates the impact of silver nanoparticles infused in blood within a stenosed artery under the influence of a magnetic field. The governing equations of continuity, mass, and heat transfer are non-dimensionalized to facilitate numerical solution. Utilizing the fourth-order Runge-Kutta method implemented in MATLAB, the non-dimensionalized equations are solved to analyze the effects of pertinent parameters on flow velocity and heat transfer in the presence of a magnetic field. The study reveals significant insights into the behavior of blood flow and heat transfer when silver nanoparticles are introduced into the bloodstream in the context of arterial stenosis and magnetic field exposure. The analysis provides valuable information on how the magnetic field and nanoparticle infusion affect flow characteristics and thermal dynamics within the stenosed artery. The rise in the volume fraction of nanoparticles slows down the nanofluid. The augmented values of magnetic parameter results in decrease in velocity but increase in temperature. The thickness of boundary layer at arterial wall decreases with enhancement in Prandtl number. The findings of this research hold promise for applications in biomedicine and medical science, offering potential strategies for enhancing therapeutic interventions in cardiovascular diseases. This study contributes to the advancement of biomedical engineering and offers avenues for the development of novel treatments and diagnostic techniques by elucidating the intricate interplay between nanoparticles, blood flow, and magnetic fields within stenosed arteries.

Pulsating magnetohydrodynamic hybrid nanofluid flow with thermal radiation, viscous dissipation, and entropy generation: Application to bio-nanomedicine

The aim of this investigation is to explore the flow phenomena in two distinct hybrid nanofluids (with micropolar fluid as the base fluid) undergoing hydromagnetic pulsation between two vertical walls, accompanied by entropy generation. This study also explores the unique heat transfer characteristics of F e 3 O 4 − Ti O 2 and A l 2 O 3 − MgO hybrid nanofluids, with the effects of Ohmic heating, thermal radiation, and viscous dissipation. The goal of using hybrid nanofluids is to improve the performance of biological systems such as nano-drug delivery in arteries, magnetic bio-separation, artificial kidneys, pressure surges, magnetofection agents, biomedical engineering, cancer treatments, and brain tumor treatment. The governing partial differential equations are transformed into a system of ordinary differential equations, which are then solved numerically by a shooting technique with a fourth-order Runge-Kutta method using a bvp4c MATLAB solver. Graphical explanations are provided for the flow variables such as the velocity, the microrotation, the temperature, the entropy generation, and the Bejan number. The heat transfer rate is an increasing function of the thermal radiation parameter; when the thermal radiation parameter increases from 1 to 3 the increases in the heat transfer rates in micropolar fluid (the base fluid), first hybrid nanofluid ( F e 3 O 4 − Ti O 2 ), and second hybrid nanofluid ( A l 2 O 3 − MgO ) are respectively 55%, 51%, and 43.86%. However, while the heat transfer rate (the Nusselt number) is a diminishing function of the Hartmann number. Further, the outcomes reveal that the hybrid nanofluids with emerging physical effects are helpful to treat cancer patients, recovering the damaged tissues and cells in arteries, and various medical-related procedures.

Mathematical modelling of HMT through porous stretching sheet using artificial neural network

The influence of heat radiation on HMT, MHD on a porous stretched sheet has been explained in this study. The governing PDEs are converted to ODEs via similarity transformations. The coefficient of skin friction, rate of HMT are calculated using the MATLAB software for various parameter values. The controlling boundary layer PDEs are turned into a system of coupled nonlinear ODEs via similarity techniques, which are numerically solved using the shooting technique in conjunction with the fourth order Runge Kutta method, and then ANN is applied to them. Validation of numerical results using ANN results. We discovered that using engineering standpoints, an ANN model may produce high-efficiency estimates for heat transfer rates. The achieved R squared values of 99% for predicting Skin friction, Nusselt number and Sherwood number coefficients highlight the remarkable effectiveness of neural network models in these predictions. This efficacy not only demonstrates the accuracy of the models but also results in a significant reduction in the computational time required compared to traditional numerical methods. Furthermore, when compared to alternative numerical approaches, the current ANN model stands out for its applicability to more complex mathematical models, because its efficiency in minimizing both time and processing capacity demands in solving such problems.

Targeted Energy Transfer Evaluation for Nonlinear Vibrations of Elastic Medium-Finite Length Beam System Using Nonlinear Energy Sink Theory

The elastic medium can usually reduce the vibration of supporting structure, and the impact of the interaction on the vibration characteristics of the structure is similar to the characteristic of nonlinear energy sink. At present, the dynamic research of beam on elastic foundation considering soil motion has been paid more and more attention. According to the modified Winkler model, the finite-depth elastic medium can be considered to the nonlinear energy sink mass, and the vibration energy dissipation capacity and parameter optimization of the elastic medium supporting finite-length beam under half sine pulse are studied. The Galerkin truncation is applied to the discretization of the governing equations. The numerical solution of the beam coupling system with simple support on the elastic medium is obtained by applying the fourth-order Runge–Kutta method. Based on this, the input energy ratio of the elastic medium dissipation is investigated. Furthermore, through the analysis and optimization of targeted energy transfer and dissipation, the dissipation effect of the finite range elastic medium on the vibration energy of its supporting beam is revealed, and the optimal parameter range of the elastic medium is proved. The results show that after adjusting the elastic medium parameters by technical means, the nonlinear energy sink can absorb most of the vibration energy of the beam quickly and effectively, and the optimal energy dissipation ratio can reach 95.16[Formula: see text]. The quantitative evaluation of the energy dissipation in elastic medium within soil–structure interaction effect is realized.

Nonlinear forced vibration of functionally graded hybrid three-phase nanocomposite toroidal shell segments reinforced by carbon nanotubes (CNTs) and graphene nanoplatelets (GPLs)

This study proposes a unique three-phase functionally graded (FG) hybrid nanocomposite material reinforcing the toroidal shell segments to investigate the nonlinear forced vibration using Reddy's higher-order shear deformation theory, von Karman geometrical nonlinearity and Stein-McElman's assumption along with the Hamilton principle. The examined toroidal shell's material consists of polymeric resin, carbon nanotubes (CNTs), and graphene platelets (GPLs) as nanocomposites. The material properties have been derived based on a modified Halpin-Tsai micromechanical model. Four distribution patterns have been studied: uniform (UD), FG-X, FG-O, and FG-V. The continuous distribution of GPLs and CNTs leads to inhomogeneous position-dependent properties throughout the shell thickness. The obtained differential equations of motion have been reduced to ordinary equations using Galerkin's technique. A multi-scales method (MSM) is used to estimate a closed-form solution representing the frequency-amplitude relation, and the state space representation is used along with numerical fourth-order Runge Kutta method (RK4) to obtain the nonlinear dynamic response of the toroidal shell. The accuracy of the current results obtained has been verified by comparing them with the relevant literature and numerical results using (RK4). In addition, the influence of both GPLs and CNTs weight fractions, nanofillers’ distribution types, longitudinal and circumferential wave numbers, elastic foundation parameters, static axial compression load, transverse excitation load, damping ratio, geometrical characteristics of the toroidal shell on the dimensionless natural frequency, nonlinear primary resonance (frequency-amplitude curve), and nonlinear dynamic response are carefully studied. The results found that combining (GPLs & CNTs) into the shell's matrix improves performance, especially in concave shells compared to convex ones. The FG-X pattern reduces peak resonance amplitude and improves natural frequency. A reduction of 12.3 % is observed for convex shells, and 21.5 % for concave shells compared to FG-O.

Heat and Mass Transfer in Unsteady Radiating MHD Flow of a Maxwell Fluid with a Porous Vertically Stretching Sheet in the Presence of Activation Energy and Thermal Diffusion Effects

The aim of this study is to analyze the effects of activation energy and thermal diffusion an unsteady MHD reactive Maxwell fluid flow past a porous stretching sheet in the presence of Brownian motion, Thermophoresis and nonlinear thermal. The non-linear partial differential equations that govern the fluid flow have been transformed into a two-point boundary value problem using similarity variables and then solved numerically by fourth order Runge–Kutta method with shooting technique. Graphical results are discussed for non-dimensional velocity, temperature and concentration profiles while numerical values of the skin friction, Nusselt number and Sherwood number are presented in tabular form for various values of parameters controlling the flow system. The present study is compared with the previous literature and found to be in good agreement.

Transient MHD Darcy-Forchheimer of Williamson-Casson flow with CCS: Application of wastewater treatment

With regards to the Cattaneo-Christoph (CCS) speculation, the thermal and mass transfer of a MHD Williamson-Casson ferrofluid flow through a permeable medium obeying the Darcy-Forchmeier law through an extended chamber or slab is carefully considered. The chemical reaction and Buongiorno nanofluid model were modified in the model to describe the nanoscale properties of liquid particles. Heat sources can be divided into two categories: linear and exponential space-dependent convection heat sources. The entire governing equations are changed to common differential terms, using a reasonable similarity change. These equations and their associated boundary conditions are calculated numerically via the fourth-order Runge-Kutta method. Comparison was made with the analytical solution in a special case, and very good agreement was reached. Economic success with deeper water purification can be attained by developing models for wastewater treatment facilities, conducting experiments and calculations on them, and making predictions for their nature in order to improve their cleaning efficiency, Prevent issues from arising during construction, operation, and reconstruction. The structural parameters of the quencher were physically modeled using mathematical modeling; It is possible to establish a cost-effective and effective filter model for industry and small settlements using the calculations and vertical filter model created in the article. Furthermore, it will be feasible to produce. Our mathematical procedure assure that: the deposit’s mass will be diminished through treated water under attractive field impacts; The magnetic field openness with nanometer-sized particles modifies the physical and synthetic properties of water particles bringing about exceptional characteristics; the water concentration will be decreased with the increase in the bending modulus; This means that in wastewater treatment the sediment mass will decrease.

Numerical Investigation of Diffusion Thermo and Thermal Diffusion on MHD Convective Flow of Williamson Nanofluid on a Stretching Surface Through a Porous Medium in the Presence Chemical Reaction and Thermal Radiation

In this paper, analyze the impact of Diffusion thermo and thermal diffusion with heat and mass transfer inherent of thermally radiant Williamson nanofluid over a stretyching surface through a porous medium under the convective boundary condition in the presence of thermal radiation and chemical reaction has been studied. The coefficients of Brownian and thermophoresis diffusions are also taken into consideration. The governing partial differential equations are reduced to a couple of nonlinear ordinary differential equations by using suitable transformation equations; these equations are then solved numerically with the use of the conventional fourth-order Runge Kutta method accompanied by the shooting technique. As a result, the effects of various physical parameters on the velocity, temperature, and nanoparticle concentration profiles as well as on the skin friction coefficient and rate of heat transfer are discussed with the aid of graphs and tables. This study has been directly applied in the pharmaceutical industry, microfluidic technology, microbial improved oil recovery, modelling oil and gas-bearing sedimentary basins, and many other fields. Further, to check the accuracy and validation of the present results, satisfactory concurrence is observed with the existing literature.

Nonlinear vibration characteristics of axially moving anisotropic membrane with gas-thermal-elastic coupling

In the actual printing production process, the printing electronic motion membrane is susceptible to transverse vibration caused by interference from air drag, temperature changes, and other external conditions, resulting in membrane wrinkles, slippage, and other phenomena in the transmission process. We studied the bifurcation and chaos movement properties of anisotropic membranes under air drag and temperature. According to D’Alembert’s theory and von Kármán’s principle, the nonlinear dynamic differential formulas of axially moving anisotropic membranes with gas-thermal-elastic coupling are established. The Galerkin method is applied to discretize the formulas to obtain the state equation of the system. Finally, numerical simulations are performed by applying the fourth-order Runge–Kutta method to analyze the bifurcation and chaos of the system in terms of orthotropic coefficient, dimensionless air drags, and dimensionless temperature. The bifurcation diagrams, Lyapunov exponent diagrams, displacement time-history diagrams, phase-trajectory plane diagrams, and Poincaré diagrams of the membrane system are obtained. The results show that the anisotropic coefficient, dimensionless air drag, and dimensionless temperature significantly impact the investigated nonlinear dynamic of the anisotropic membrane, which provides a theoretical basis for production efficiency and high-quality printing equipment.

Effect of Stratified Flow on the Vibration of Anchor Cables in a Submerged Floating Tunnel

This study investigates the vertical-type submerged floating tunnel with anchor cables. Based on the characteristics of the anchor cables, the anchor cables are simplified as a nonlinear beam model with hinged ends. Disregarding the axial displacement of the tunnel body, the loads will cause displacements in the x and z directions of the tunnel body. The vibrations of the anchor cables are decomposed into three directions, and the parameter excitation at the connection point between the anchor cables and the tunnel body is taken into account. The equations of motion for the three degrees of freedom of the anchor cables are established using Hamilton’s principle, and then the three equations are solved using the Galerkin method and the fourth-order Runge–Kutta method. The basic characteristics of an internal wave stratified flow acting on the anchor cables are considered, as well as the influence of the incident angle of the ocean currents on the three degrees of freedom of the anchor cables. The results indicate that (1) stratified flow weakens the first- and third-order vortex-induced vibrations of the anchor cables while enhancing the second-order vortex-induced vibrations. When considering the parameter excitation of the anchor cables, the first- and third-order vibrations are weakened, while the second-order vibration remains significant; (2) the first-order vibration of the anchor cables reaches its maximum value when the transverse oscillation frequency of the tunnel body is twice its natural frequency, and the second-order vibration of the anchor cables reaches its maximum value when the transverse oscillation frequency of the tunnel body is twice its natural frequency; (3) the downstream vibration of the anchor cables increases with the increase in the incident angle of the ocean currents, the cross-flow vibration of the anchor cables decreases with the increase in the incident angle of the ocean currents, and the axial vibration of the anchor cables reaches its maximum value when the incident angle of the ocean currents is 60 degrees; (4) stratified flow weakens the lock-in phenomenon of the anchor cables, and the influence of the 1/2 stratified flow on the vibrations of the anchor cables is greater than the influence of the 1/2 stratified flow.

High-efficiency and low-hazard artillery recoil reduction technology based on barrel gas reflection

Reducing recoil as well as reducing muzzle hazards are important issues in artillery design. This paper presents a barrel gas reflection method for the artillery aiming for efficient recoil reduction while reducing muzzle hazards. The launching process is modeled by coupling the interior ballistic equations and the flow equations of the barrel gas reflection device. The fourth-order Runge–Kutta method was used to solve the model, and the LHS method as well as the Kriging model was used to establish a mapping relationship between the parameters and the effect. To validate the proposed model, shooting experiments are carried out on a 30 mm caliber artillery. The maximum error between the experiment and simulation results was 5.32%. The experiment has demonstrated that the barrel gas reflection method increases the recoil reduction efficiency of artillery by 44.54% and reduces the muzzle hazard by 52.18%. Finally, the barrel gas reflection method can produce effects with the muzzle device at the same time, and it has little effect on the velocity of the projectile muzzle velocity, and it provides a new way of thinking for the development of future artillery recoil reduction technology.

New Analytical and Numerical Solutions for Squeezing Flow between Parallel Plates under Slip

In this article, the effects of physical flow parameters on squeezed fluid between parallel plates are explored through the Darcy porous channel when fluid is moving as a result of the upper plate being squeezed towards the stretchable lower plate, such as velocity slip, thermal slip, solutal slip, thermal stratification parameter, solutal stratification parameter, squeezing number, Darcy number, Prandtl number, and Schmidt number. The governing equations are transformed into a nonlinear ordinary differential equation using the appropriate similarity transformations. The resulting equations are solved by using the perturbation iteration method (PIT) to produce a convergent analytical solution with high accuracy. The phenomena of the squeezing fluid as the plates are moving apart and when they are coming together are illustrated using the resulting analytical solutions. Plots are used to discuss the significant effects of physical parameters on velocity, temperature, and fluid concentration profiles. The skin friction coefficient and Nusselt Sherwood values have graphical interpretations that are listed. For strong velocity slip parameters, the results demonstrate the existence of a minimum velocity profile close to the plate and a growing velocity profile distant from the plate. Additionally, as the slip effects rise, the fluid temperature and concentration both considerably drop. The results of the fourth-order Runge-Kutta method (RK4M) and the presented analytical solutions provided are in excellent agreement.